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Abstract 

Systematic variations of average observational characteristics and correlation properties 
■ between different parameters of gamma-ray bursts (GRBs) with time from 1991 April to 

1994 September are revealed. It is hard to explain the observed long-term variability 

by variations of experimental conditions. The variability of GRB properties with the time 
t^j- ' scale of months to years, together with the similarity between GRBs, solar hard X-ray flares 

and terrestrial gamma-ray flashes, may indicate an origin of GRBs, at least partly, within 
^] \ the heliosphere. Large-voltage and high-temperature pinch plasma columns produced by 

£2 ' disruptive electrical discharges in the outer heliosphere can generate hard X-ray and 7-ray 

flashes with energy spectra and spectral evolution characters consistent with that observed 

in GRBs. 
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1. Introduction 

^ ■ GRBs are intense, brief flashes of non-thermal hard X-rays and 7-rays from random directions and at 

unpredictable times. GRBs have been extensively studied since their discovery (Klebesadel, Strong & 
Olson 1973). Several satellites were launched to observe the bursts and a large number and wide variety 
of models were put forwards to explain their origin. The Burst and Transient Source Experiment 
(BATSE) on board the Compton Gamma Ray Observatory (CGRO) was launched in April 1991. 
BATSE, with an unprecedented sensitivity, ~ 1 x 10 -7 erg cm -2 , has revolutionized the study about 
the nature of these objects and raised more debates on it as well. 

After more than twenty five years of the discovery of GRBs, their origin still remains a continuing 
puzzle. The most significant characteristic of these transient events is their complexity in morphol- 
ogy. The duration of GRB ranges from a few milliseconds to hundred of seconds and the temporal 
structure displays complicated patterns, of which basic features are still difficult to be summarized. 
In all timescales, some bursts have complex profiles, rich in fluctuating structures and some smooth, 
sometimes single structure profiles. None of the models suggested so far for GRBs can provide a 
satisfactory explanation to their main features in temporal and spectral characteristics. 

The observed bursts are distributed isotropically on the sky, that exclude distant galactic disk 
populations from the GRB sources. The distributions of GRB intensities (peak counts C or peak 
flux P), integral IgN — lgC or lgA r — lgP diagrams, are consistent with the —3/2 slope (expected 
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for a homogeneous sample of monoluminosity in Euclidean space) for bright bursts, but considerably 
absent for bursts at low intensities. The incompatibility of the intensity distribution with a uniform 
distribution throughout space rules out local galactic disk sources. The only place left for GRBs 
at the Galaxy is at an extended galactic halo with distances d ~ 50 — 100 kpc (Hartmannm 1994). 
The intensity distribution is widely accepted as a strong support to an extragalactic origin that the 
deviation from the —3/2 power law is due to cosmological effects (Mao & Paczyhski 1992; Piran 1992; 
Dermer 1992). The debate on the distance scale to GRBs is now mainly between the galactic halo 
and cosmological distances (eg. Lamb 1995; Paczyhski 1995). 

Some average properties of GRBs detected by BATSE, e.g. duration and spectral hardness, show 
systematic variations with time from 1991 April to 1994 September (Li 1996). These variations are 
difficult to be explained with experimental effects. The analysis procedures and more results on the 
variability of properties of correlation between different GRB parameters are presented in section 2 of 
this paper in detail. The existence of GRB long-term variability with a time scale of months to years 
indicates a heliospheric origin. In section 3 a model of GRBs originating from disruptive discharges in 
instable plasma of the outer heliosphere is described. A study of two kinds of 7-ray flashes, solar flares 
and terrestrial 7-ray flashes, which were also detected by BATSE and certainly originated within the 
heliosphere, is made in comparison with GRBs in section 4. 

2. Long-term variability 

BATSE, having been continually detecting GRBs at a rate of about one event per day since April 
1991, makes it possible to study long-term variability of GRB observational properties. The third 
BATSE (3B) Catalog of gamma-ray burst (Meegan et al. 1996) consists of 1122 7-ray burst triggers 
and covered the time interval from 1991 April 21 until 1994 September 19. To check through the 3B 
catalog to see if there exists systematic change with time in values of a burst parameter, we divided 
the observation time (TJD from 8367 to 9614, where the truncated Julian Date TJD=JD-2440000.5) 
into a few epoch periods with the same length and calculated the average value of the parameter for 
each period. 

2.1 Duration and hardness 

The parameter Tqq is used to measure the duration of a burst. T90 is the time interval in which the 
integrated counts from the burst increases from 5% to 95% of the total observed counts. Firstly we 
calculated the average T and standard deviation ot of all the values of T90 listed in the 3B catalog, 
and then divided the observation time into four periods and for each of the four periods calculated 
the average, < T90 >, from all the n bursts in each period for which T90 are available. The result is 
shown in Fig. 1(a). The statistical error of < T90 > is estimated by o^l yfn. One can see from Fig. 1(a) 
that < T90 > evidently decreases with time, from 41.4 ± 3.6 s during 91.4 — 92.2 to only 25.0 ± 3.4 s 
during 93.11 — 94.9, at a rate , relative to T, of (—18.3 ± 5.6)% yr -1 , which was derived by a weighted 
linear least-squares fit with a reduced \ 2 = 0-21 (the reduced \ 2 ls the fitting x 2 over the degrees of 
freedom). The linear regression has 3.3<r significance for the null hypothesis that T90 and observation 
time are uncorrelated. 

The variability of T90 has also been studied by the correlation analysis without binning. The 
correlation coefficient for a sample of n pairs of x, y values is defined as 




i=i 
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Figure 1: Average burst duration < Tg > against observation time of 3B bursts, (a) For all the 
bursts for which T90 are available in 3B catalog; (b) For bursts in 3B catalog with trigger thresholds 
Cmin < 75 cts on the 64 ms timescale. 



where v{i) = (x(i) — x)(y(i) — y). For the total n = 834 bursts in 3B catalog for which T90 are 
available, the value of the correlation coefficient between T90 and the observation date t was calculated 
as r(T$o,t) = —0.12. The significance of a value of r to test the hypothesis that T90 and t are 
uncorrelated can be evaluated by the statistic u = ry/Jn — 2)/(l — r 2 ) which satisfies a Student's t- 
distribution with n — 2 degrees of freedom under the null hypothesis (Siegel 1956). We found u = —3.4. 
As n = 834 3> 1, the statistic u here follows the standard normal distribution, u = —3.4 indicates a 
result with 3.4<r significance, which is consistent with the result from regression analysis (3.3<r) and 
the probability that the observed variability of T90 comes from the statistical fluctuation is less than 
3.2 x 10~ 4 . 

In the 3B catalog four energy fluences Fi,F2,Fs and F4, in units of erg cm -2 and covering the 
energy ranges 20-50, 50-100, 100-300 keV and E > 300 keV respectively, are given for 867 bursts. 
The average fluences of bursts in each of four periods were calculated for the four energy channels 
separately and shown in Fig. 2 (a). The fitted change slopes from Fig. 2 (a) are listed in Table 1. One 
can see from Fig. 2(a) and Table 1 a steady trend that the average fluence increases quicker for higher 
energy channel. 



Table 1: Change rates of average fluences 



Energy ran£ 


re (keV) 


20 - 50 


50 - 100 


100 - 300 


> 300 


Change rate 


(% jO 


10.2 ± 10.2 


6.2 ±7.3 


17.8 ±9.7 


29.0±15.8 



The ratio H24 = F2/F4 is used to measure spectral hardness. In statistics of H24 a few events 
with very small F4 seriously influence the average and disturb the general trend of variability much. 
Only three isolated bursts in total 690 bursts with F4 > have H24 > 40, we picked them out and 
calculated < H24 > and its standard deviation for residual 687 events for each period. The results 
shown in Fig.2(b) reveal that the energy spectra of GRBs in 3B catalog became harder with time 
(< H 2 4 > decreased, from 1.15 ± 0.17 during 91.4 - 92.2 to only 0.44 ± 0.18 during 93.11 - 94.9, at a 
rate of (-35.1 ± 11.9)% yr" 1 ). 
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Figure 2: Evolution of burst spectral property, (a) Average energy fluences against observation time 
for 3B bursts. Fi,F2,F$ and F4 are fluences covering the energy ranges 20-50, 50-100, 100-300 keV 
and E > 300 keV respectively; (b) Average spectral hardness ratio against observation time for 3B 
bursts. The hardness ratio H24 = F2/F4, from all the bursts for which F2 and F4 are available and 
F2/F4 < 40. 

2.2 Experimental effects 

Since March of 1992 there are gaps in the data of GRB intensity, fluence or duration due to the errors 
of CGRO tape recorders, the problem has significantly reduced since March of 1993 (Meegan et al. 
1996). The number n of events for which T90 can be measured and that of total detected GRBs, N, 
for each period, are listed in Table 2. A less value of the ratio n/N reflects greater effect of data gaps. 
As the change behaviour of < T90 > and < H24 >, shown in Fig.l(a) and Fig.2(b), is quite different 
with that of n/N, shown in Fig.3(a), it is hard to interpret the systematic variations of duration and 
hardness ratio by the effect of data gaps. 

The observed variation of burst duration can not be explained by the variation of background 
either. The average of trigger thresholds C m i n on the 64 ms timescale, set by command to a 5.5cr 
deviation above background, was calculated for each period and also listed in Table 2 and shown in 
Fig.3(b). The values of < C m j n > have not shown a systematic change with time. 



Table 2: Statistics of 3B bursts in four periods 



Period 


Number of Bursts 




<T 90 > 


< H 24 , > 


(yy.mm.dd) 


n 


N 


n/N 


(cts/64ms) 


(s) 




91.04.21-92.02.27 


216 


249 


0.87 


69.9±1.5 


41.4±3.6 


1.15±0.17 


92.02.28-93.01.04 


168 


271 


0.62 


70.6±3.5 


35.6±4.1 


0.80±0.21 


93.01.05-93.11.12 


208 


297 


0.70 


66.7±0.9 


33.2±3.7 


0.66±0.19 


93.11.13-94.09.19 


242 


305 


0.79 


72.7±2.6 


25.0±3.4 


0.44±0.18 



For further inspecting the effect of background variation on < T90 >, the C m in distribution for 
bursts in the 3B catalog was made and shown in Fig. 4(a). Most values of C m i n are smaller than 75 
cts/64ms. The bursts were grouped into four subsets with C m i n < 65, 65 < C m i n < 68, 68 < C m i n < 75 
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Figure 3: Variations of BATSE data gaps and backgrounds, (a) Efficiency ratio n/N for measurement 
of Tgo against observation time, where N is the total number of detected bursts during a period and 
n the number of bursts for which Tgo can be measured during the same period, (b) Average trigger 
threshold on the 64 ms timescale against observation time. 
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Figure 4: Distribution of trigger thresholds on the 64 ms time scale of 3B catalog; (b) Correlation 
between duration and trigger threshold of 3B bursts. 



and C m in > 75 (in units of cts/64ms) respectively. The average values of Tgo for each subsets were 
calculated and shown in Fig.4(b). There is no visible dependence between < Tgo > and C m i n except 
such detected bursts under highest trigger thresholds of C m i n > 75 cts/64ms. Fig. 1(b) shows that the 
average duration < Tgo > for bursts having C m i n < 75 cts/64ms also decreases with time at a rate of 
(—17.1 ± 5.2)% yr _1 , consistent with that derived from all bursts. Therefore, the observed duration 
variability can not be interpreted with background variation. 

Another experimental condition with secular change is spacecraft altitude. Fig.5 shows CGRO 
altitude history. After comparing Fig.5 with Fig.l and Fig. 2, one can see that the observed long-term 
variations of duration and hardness can not be due to the altitude variation either. 

There exists a positive dependence of burst fluences against their durations. Fig.6(a) shows the 
correlation of lgT and lg < Tgo > in 3B events where the regression line has a slope of A = 0.52 ±0.02. 
If the decrease of average duration is caused by decrease of detection efficiency for longer events by 
some unknown reasons, a decrease of average fluence with time should be expected. But the measured 
results for fluences show an opposite trend, see Fig.6(b), indicating that the detected changes should 
be intrinsic. 
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Figure 5: CGRO altitude vs. time. 




Figure 6: (a) Correlation between fluence and duration of 3B bursts; (b) Average fluence against ob- 
servation time, solid lines represent the measured results and dashed line is the expectation calculated 
from the regression line for < Tgo > in Fig. 1(a) and correlation between F and T90 and normalized at 
TJD=8367. 
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2.3 Correlation properties 

The relationship between two parameters of GRBs reflects the nature of GRB production process. 
Long-term variability also found for different correlation property of GRBs in 3B catalog, which is 
more difficult to be interpreted with experimental effects (Li 1997). There are total 807 bursts in the 
3B catalog for which both F and Tgo are available. Calculating the correlation coefficient between 
lgT and lgT 90 for the 807 bursts we get r(lgT,lgT 90 ) = 0.67 ± 0.02. The statistical error of the 
correlation coefficient is estimated by the bootstrap technique (Efron 1979; Diaconis & Efron 1983). 
To check through the 3B catalog to see if there exists systematic change with time in the relationship 
between F and Tgo, we divide the observation time t (TJD from 8367 to 9614 for 3B catalog) into 
five epoch periods, letting each period has nearly equal number of bursts for which both F and 
T90 are available, and calculate the correlation coefficient r(\gF, lgTgo) and its statistical error for 
each period. Fig. 7(a) shows the result. From Fig. 7(a) one can see that the correlation coefficient 
r(\gF, lgT 90 ) increased from 0.60 ± 0.04 during 91.111 - 91.329 to 0.71 ± 0.03 during 94.154 - 94.262 
at a rate A(r,t) = 0.038 ± 0.017 yr" 1 , which was derived by a weighted linear least-squares fit. The 
linear increasing of correlation coefficient r(lgT, lgTgo) with observation time t has 2.2<r significance 
for the null hypothesis that r(lgT, lgTgo) and observation time t are uncorrelated. 

Now we study the variability of correlation coefficient between spectral hardness and duration. The 
ratio H32 = F3/F2 is used to measure spectral hardness. As events with very small F2 considerably 
disturb the general trend of variability of H32, we pick a couple of events with H32 > 10 out in the 
statistics through this work. There is a general trend of hardness with short events being harder (Deza- 
lay et al. 1992; Kouveliotou et al. 1993a), but the variation behaviour of hardness against duration 
for short bursts is quite different with that for long bursts: the correlation coefficient r(H^2, lgTgo) 
between H32 and lgTgo for 583 long bursts with Tgo > 2 s in the 3B catalog is 0.01 and that for 178 
short bursts with Tgo < 2 s is —0.11. No visible change of the correlation coefficient between H32 and 
lgTgo with time has been found for long bursts, the slope of the regression line of r (H32, lgTgo) on t 
fitted to the values of five epochs being 0.005 ± 0.039 yr _1 , whereas evident variation existed in short 
bursts, shown in Fig. 7(b). The correlation coefficient r(H^2, lgTgo) f° r bursts with Tgo < 2 s increased 
from -0.52 ± 0.11 during 91.111 - 91.324 to 0.16 ± 0.13 during 94.43 - 94.262. The regression line in 
Fig.7(b) has a slope of 0.239 ± 0.057 yr -1 , showing a 4.2<r long-term variability. 

2.4 Bimodality 

The bursts duration distribution can be divided to two sub-groups according to T 9 o: Long bursts with 
Tgo > 2 s and short bursts with Tgo < 2 s (Kouveliotou et al. 1993b; Meegan et al. 1996). We 
found that the variability behaviour of observed GRB parameters is also bimodal: duration becoming 
shorter and spectrum harder is only observed for long bursts in the 3B catalog, but not for short ones. 
On the other hand, the variability of correlation properties between different parameters, e.g. F and 
Tgo, or H32 and Tgo, is revealed in short bursts more obviously than that in long bursts. Fig. 8 show 
time variation behaviour of the parameters (duration, hardness ratio and fluence), and the correlation 
coefficients between fluence and duration and between hardness and duration for long bursts and for 
short bursts in the 3B bursts separately. The change slopes derived by the linear fitting procedure are 
listed in Table 4. For short bursts both r(lgT, lgTgo) an d r(Hs2, lgTgo) increased with time evidently 
with Act significance, but for long bursts the change rate of r(lgT, lgTgo) was only one third of that 
for short bursts and no visible change found for r(Hs2, lgTgo). 
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Figure 7: Evolution of correlation coefficients, (a) Correlation coefficient between IgF and lgTgo 
against observation time for 3B bursts, (b) Correlation coefficient between hardness ration H32 and 
lgT 90 against observation time for short 3B bursts with T go < 2 s. 



Table 3: Change slopes of GRB properties 



Property 


Change 


slope (yr i ) 




short bursts (Tgo < 2 


s) long bursts (Tgo > 2 s) 


Tgo 


(-2.5 ±5.8)% 


(-15.7 ±5.0)% 




(12.0 ±27.2)% 


(-32.3 ± 12.9)% 


F 


(16.3 ± 14.6)% 


(19.7 ± 13.2)% 


r(lgF,lgT 90 ) 


0.20 ±0.05 


0.064 ± 0.028 


r(# 32 ,lgT 9 o) 


0.24 ± 0.06 


0.005 ± 0.039 



2.5 Short bursts 

From Fig.8(i) and (j) and Table 3 one can see that long-term variability of the correlation coefficient 
r(Hz2, lg Tgo) only show in short bursts with Tgo < 2 s, none in long bursts. Now we further inspect 
the correlation properties of short bursts. 

There are 745 events in the 3B catalog for which both T90 and F are available. To show the 
relationship between hardness and duration, we divide the range — 200 s of T90 into 22 bins, letting 
each bin includes nearly equal number of events, and calculate the average of Tgo, average of H32 and 
their standard deviations for each bin. Fig.9 shows the result. It is noticeable in Fig.9 that a turn in 
the correlation behaviour between H32 and Tgo occurs at Tgo c± 0.55 s as well as at 2 s. From Figure 10 
(a) and (b), which show the H32 — lgTgo distribution of 373 bursts during observation epoch t between 
8367 and 9037 TJD and that between 9038 and 9614 TJD separately, one can see that the correlation 
behaviour between H32 and Tgo of the brief bursts (Tgo < 0.55s) is evidently changed with time. 

Similar phenomena can also be found in the correlations between H41 and T90 and that between the 
peak flux Pg4i the maximum flux in 50 — 300 keV (ph. cm^s" 1 ) integrated over 64 ms timescale, and 
Tgo for the short bursts. The correlation coefficient r(H^i, lgTgo) evidently increased with observation 
time at a rate A = 0.36 ± 0.08 yr" 1 for the brief bursts shown in Fig.ll(a), but no visible variation is 
found in Fig.ll(b) for the medium bursts (A = 0.008 ± 0.12 yr" 1 ). For the brief bursts the correlation 
coefficient of v(Pq4, lgTgo) is positively correlated with the observation time (the slope of regression 
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Figure 8: Evolution of average values of parameters Tqq,H24,,F and correlation coefficients 
r(lgF, lgTgo), r(i?32, lgTgo) for two types of GRBs in the 3B catalog. (a),(c),(e),(g),(i) for bursts 
with T 90 < 2 s, (b),(d),(f),(h),(j) for bursts with T 90 > 2 s. H 24 = F 2 /F 4 and H 32 = F 3 /F 2 . A couple 
of bursts with H 2 4 > 40 were excluded from the statistics in (c) and (d), and that with H 32 > 10 in 
(i) and (j). 
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Figure 9: Distribution of H32 — lgTgo for all the events in 3B catalog during the period between 8367 
and 9614 TJD. 
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Figure 10: H32 — lgTgo distributions, (a) For bursts during the period between 8367 and 9037 TJD. 
(b) For bursts during the period between 9038 and 9614 TJD. 
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line A = 0.26 ± 0.06 yr^ 1 in Fig.ll(c)), but for the medium bursts a negative correlation between 
v(Pq4, lgTgo) and the observation time appears with A = —0.15 ± 0.10 yr _1 shown in Fig. 11(d). 

To get a free of binning estimation of the significance level for the time variability of the correlation 
between x and y from a sample of n points (x,y,t), the equation of the line of regression of x on 
t, x = A x + B x t, and that of y on t, y = A y + B y t, are firstly derived by the method of least 
squares. For each burst i (i = 1, ...,n) with parameters (x(i),y(i),t(i)) the following quantity v{i) = 
(x(i) — xo(i))(y(i) — yo(i)) with xq(i) = A x + B x t(i) and yo(i) = A y + B y t{i) is calculated. The 
correlation coefficient between v and t, r(u,t), for the sample of n pairs of v,t- values now can be 
calculated, the standard deviation of the derived r(u, t) can be evaluated by the bootstrap technique 
and the relative portion £ of the correlation coefficient less or equal to (or greater or equal to if 
r(u, t) < 0) among the bootstrap samples can be taken as an estimate of the significance level. The 
obtained values of r(v,t) of the parameter pairs (#32, lgTgo), (H41, lgTgo), and (P@4, lgTgo) are listed 
in Table 4 for the two types of burst separately, where the standard deviations and significance levels 
are derived from 10 7 bootstrap samples for each case. 



Table 4: Significance levels £ of correlation variations 



Correlated Parameters 


T 90 < 0.55s 


0.55s< T 90 < 2s 


(x, y) 


r(u,t) £ 




H32, lgTgo 


0.30 ±0.09 2.1 x 10" 4 


-0.03 ±0.10 0.38 


#41, lg?90 


0.39 ±0.07 9.3 x 10~ 5 


-0.07 ±0.11 0.26 


^64, lgTgo 


0.36 ±0.07 4.8 xlO" 6 


-0.16 ±0.09 4.1 x 10" 2 



One can see that the three correlation coefficients listed in Table 4 have quite different variation 
behaviour between the two types of burst and significant (near to or great than 4a) variations appear 
in the correlation properties of the brief bursts (Tqq < 0.55 s). 



3. A discharge model 

Alfven (1981) has stressed the importance of studying electric current circuits in plasmas to under- 
standing phenomena occurred in the magnetosphere till galactic dimensions, and held that many of 
the explosive events observed in cosmic physics are produced by exploding electric double layers. In 
a circuit with inductance L and current I in plasma, electric double layers can be produced, in which 
large potential drops may be built up over distances of the order of some tens of the Debye lengths. 
A disruption of the current in a layer by plasma instability will cause an explosion: the layer voltage, 
V = L(dl/dt), may exceed the normal value by several orders of magnitude and the magnetic energy 
stored in the circuit being suddenly released in the layer, will cause the plasma column to collapse 
to a much smaller radius by inward magnetic pressure of the discharging current (Z pinch) (Krall & 
Trivelpiece 1973) and the resulting compressed, large-voltage and high-temperature discharge column 
will radiate energetic photons from the inverse Compton scattering between the high-energy electrons 
and thermal photons. In an isotropic thermal emission with mean photon energy hv, the mean energy 
of a scattered photon from an electron with an energy e is E = ^{^i) 2 hv (Ginzburg & Syrovatskii 
1964). Let the temperature of a pinch discharge column is T (K), the differential photon density of 
the thermal emission at frequency v is n u (T) = B U (T)[1 — exp(— t(v))], where B U (T) is the photon 
density of a black body and the optical depth 



t{u) = 4.1 x 10- 23 T- 3 - 5 (-pr 3 (l - e~w)g f N 
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Figure 11: Correlation coefficient vs. observation time. For each figure the observation time period 
of 3B catalog is divided into five epoch periods by letting each period has nearly equal number of 
available events. The correlation coefficient r is calculated for each bin and its standard deviation is 
estimated by the bootstrap technique. A couple of isolated events with relevant parameter > 150 
or Pq4 > 15 are excluded from statistics, (a) r(i?4i, lgTgo) vs. observation epoch for Tgo < 0.55s. (b) 
r(i?4i, lgTgo) vs. observation epoch for 0.55s< Tgo < 2s. (c) r(P^,\gT§o) vs. observation epoch for 
Tgo < 0.55s. (d) r(P64,lgTgo) vs. observation epoch for 0.55s< Tgo < 2s. 
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with the emission measure / N 2 dl being in cm 5 and average Gaunt factor g taken as the approxima- 
tion (Kulsrud 1954) 

_ _ V3 y/0.76 + hv/kT + 0.87 
9 ~~ n ^0.76 + hiTJkT - 0.87 

For thermal emissions of temperature T = 10 6 K and optical depth To = 0.1 at hv = lkT, the mean 
photon energy ~ 72 eV, then the mean energy of a scattered photon from an electron with e = 100 
MeV can be calculated as ~ 3.7 x 10 3 keV. Therefore, the large voltage and high temperature discharge 
column is an effective radiator of hard X-rays and soft 7-rays. 

The existence of any intrinsic variability of GRB properties with a time scale of months to years 
can not be interpreted with either the galactic or cosmological origin. Only the decrease in solar 
activity during the 3B observation period remains to be a reason, and thus GRBs, or at least a 
part, should be originated within the effective region of the solar wind. GRBs may be possible from 
exploding electric double layers in instable plasma within the heliosphere (Li 1996). With respect 
to their isotropic spatial distribution, GRBs should be produced far away from the Sun. A distance 
of about 100 AU can satisfy the constraint from observed GRB isotropy distribution (the measured 
dipole moment cos 9 ~ 0.01) and that from triangulation and direct positioning for a few bursts. At 
the solar-wind terminal shock, near the heliospheric boundary and surrounding the solar system, lies 
an abrupt discontinuity, turbulent and complex structure (Suess 1990), at which various intriguing 
phenomena may occur, GRB might be one of them. The typical energy released in a GRB produced 
at a distance of 100 AU is about 10 26 / ergs, where / is the beaming solid angle fraction. If emission 
is beamed into an angular region of ~ 1 deg 2 , the typical energy of GRBs will be only ~ 2 x 10 21 ergs, 
equivalent just an average thunderstorm (Battan 1961). 

A burst with total released energy 2 x 10 21 erg emits ~ 2 x 10 28 photons (60 keV/ph assumed). 
If each discharge electron produces one burst photon through the inverse Compton scattering with 
thermal photons, a discharging current of ~ 2 x 10 s A can produce the burst emission with a 10 s 
duration. A current flowing in the heliospheric plasma may contract by the magnetic confinement 
and form a plasma cable with much larger density than the surroundings (Alfven 1981). In the 
disruptive discharge the cable is further pinched into a very narrow column. The interaction length 
of the Compton scattering between energetic electrons and thermal photons in a plasma column with 
diameter d=lm and temperature 10 6 K, which is pinched from a initial current cable with diameter 
do = 0.1 AU= 1.5 x 10 12 cm and density Nq = 1 cm~ 3 , is only ~ 10 km. The density of the pinch 
column can be estimated as A ~ 2 x 10 20 cm" 3 and optical depth To ~ 0.1. 

The GRB spectra accumulated by the Spectroscopy Detectors (SD) of BATSE are described well 
by the following empirical Band et al. (1993) model 

N(E) = Ae a exp(-E/E ) (a-(3)E >E 
= A[(a - /3) J E ] Q_/3 exp(/3 - a)E p (a - (3)E < E 

with a ~ — ~ —2. The peak power energy E p k = (2 + ol)Eq. Our calculations show that 
7-ray flashes produced by large voltage and high temperature discharge plasma columns can also 
be described well by this model (Li & Wu 1997). We assumed a simple discharge model that an 
initial discharge voltage Vq is produced at time t = 0, the discharging current firstly has a Gaussian 
rise, i(t) oc exp(— (^^) 2 ) till to its maximum at t = t m and then falls exponentially with a decay 
constant a,i(t) oc exp(— *~^ m ); the potential difference falls proportionally to the total quantity of 
charge transferred by the current, V(t) = Vq — B jQi(t)dt. At time t, the kinematic energy of a 
discharge electron increases lineally along the path from at the cathode to V(t) at the anode, the 
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Figure 12: An example of discharge model, (a) Time histories of discharge voltage (solid line) and 
current (dashed line), the voltage is in units of 10 6 V, current in arbitrary units, time is scaled by 
the exponential decay constant of the current, (b) The calculated time-averaged spectrum (crosses) 
from a discharge column with temperature T = 1 x 10 6 K, optical depth at lkT tq = 0.1 and time 
histories of voltage and current shown in (a), observed over a scattering angle interval 0°- 0.5°. The 
flux is in arbitrary scale. The solid line shows the Band et al. spectrum with parameters a = —1.12, 
(3 = -2.28 and E pk = 826.4 keV. 

observed radiation is a summation of the scattered photons from electrons with kinematic energies 
between and V(t). For example, let Vq = 100 MV, a± = 0.4a, total discharge time D = 2.5a, the 
voltage at the end of discharge Vb = 5 MV, the corresponding time histories of discharge voltage and 
current are shown in Fig. 12(a), where the time is scaled by the current decay constant a. Differential 
scattered photon spectra from electrons with different energies and isotropic thermal emission with 
T = 10 6 K and tq = 0.1 were calculated by the Klein-Nishina formula. Because the plasma instability, 
the discharge column may be twisting, flickering and branching, an observed spectrum should be an 
average over a certain interval of scattering angle. As an example, the crosses in Fig. 12(b) show the 
expected spectrum observed by SD of BATSE for scattered photons with scattering angle < 0.5°, 
the error bars were estimated by assuming the total number of recorded photons of E > 20 keV 
being 2 x 10 4 , a E~ 2 background spectrum with a flux of 1 x 10~ 2 cm s _1 keV _1 at 100 keV and 
in considering the SD full-energy peak efficiency (Fishman et al. 1989). Fitting the Band et al. 
model to the simulated spectrum in the region of E > 20 keV we got the solid line in Fig. 12(b) with 
a = —1.1, (3 = —3.2, E p k = 748.0 keV and the reduced \ 2 = 0.59, showing the discharge model can 
produce a typical GRB spectrum. 

Producing Band et al. spectra is a general property of the discharge mechanism. We calculated 
100 spectra with different values of the discharge model parameters. In our calculations, the initial 
voltage Vo was sampled uniformly from the region of 50 — 10000 MV; the temperature of thermal 
emission, T, sampled from 5 x 10 5 — 5 x 10 8 K and optical depth at lkT, tq, from 0.01 — 1.0 uniformly 
on a logarithmic scale; the current decay constant was taken as unit, a = 1; the current rising constant 
a\ was sampled from 0.01 — 0.7; the voltage at the end of a discharge sampled from 0.5Vo — 5 MV; 
the center of scattering angle region considered was taken randomly from 0°- 0.5°and the width of 
the interval from 0.1°- 0.5°. For each parameter set, we calculated the expected spectrum from the 
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Figure 13: Distribution of a and from the fit of Band et al. model to calculated spectra for 100 
samples of discharge model. 

discharge process and fitted the Band et al. model to. Most spectra fit well with the model, the 
average reduced \ 2 = 1-2 ±0.3. Fig.13 is a plot of a vs. E p k from our calculations for 100 samples of 
discharge model, which shows a distribution quite similar with that from observed GRBs (see Preece et 
al. 1996, Fig. 6). More calculations show that a summation of several spectra from different discharge 
columns with different parameters can also be well described by the Band et al. model. 

A hard-to-soft spectral evolution has been observed for many bright GRBs (Ford et al. 1995; 
Briggs 1995): the peak energy, E p k, softens over the burst. Liang & Kargatis (1996) further found 
that for a sample of bursts consisting of well resolved, isolated pulses E p k decreases exponentially 
with photon fluence (running time integral of the flux), and that for many multi-pulse bursts, the 
exponential decay constant is invariant from pulse to pulse. The above properties can be interpreted 
easily in the frame of discharge mechanism. For example, we calculated time-resolved spectra for the 
discharge model described by Fig. 12. To focus on the evolution of E p k, we followed Liang h, Kargatis 
(1996) to fix the two slopes throughout the burst by fitting the Band et al. model to the time- 
average spectrum Fig. 12(b), the obtained time history of E p k, Fig.l4(a), shows a typical hard-to-soft 
evolution. In Fig. 14(b) E p k is plotted against photon fluence (the relative errors of E p k are simply 
taken as 10%, a typical value for bright GRBs of BATSE), showing an exponential decay. It is found 
from calculations that the decay constant of E p k versus fluence is mainly determined by the velocity 
of discharge voltage decaying, i.e. by the parameter B of discharge circuit. Therefore, the similarity 
of the exponential decay constants of the subsequent pulses to that of the leading pulse in samples of 
multi-pulse bursts can be explained as such pulses are from repeat discharges of the same circuit for 
the leading pulse. The softening trends in E p k is not a universal property for the discharge model. 
For the case of both the initial discharge voltage and temperature being high, the cross section of 
Compton scattering falling with the increase of electron energy (Akhiezer & Berestetskii 1965) causes 
Epk has not considerable change or even increase with time. Ford et al. (1995) have already observed 
such evolution trend in a few bursts. More dedicated model of discharge may explain the phenomenon 
of envelope decaying observed in long multi-pulse bursts (Ford et al. 1995). 
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Figure 14: Spectral evolution expected for the discharge model described by Fig. 12. (a) Time histories 
of E p k (keV) (crosses) and photon flux (dashed line, in arbitrary units). The time is scaled by the 
decay time constant of discharging current, (b) E p k versus photon fluence in units of total fluence. 
Relative errors of E p k are taken as 10%. Dashed line shows an exponential decay. 

Another natural deduction of the discharge model is the durations becoming shorter with increasing 
energy. For the discharge column described by Fig. 12 we calculated the light curves in four energy 
channels: 20 — 50, 50 — 100, 100 — 300 and > 300 keV, and their autocorrelation functions respectively. 
The calculated autocorrelation function decreases more rapidly with time lag in the higher energy 
channels compared to the lowest one, as shown in Fig. 15, which is consistent with that observed 
for GRBs (Link et al. 1993). We also calculated the cross-correlation functions between channel 1 
(20 — 50 keV) and 3 (100 — 300 keV) and found a maximum at the time lag of —0.01 (in units of the 
time constant a of current decay), indicating a time delay between soft and hard emissions, which is 
also observed in GRBs (Chipman 1993; Cheng et al. 1995). 

4. Discussion 

The morphological diversity, the most striking feature of GRBs, is not a difficulty for the discharge 
mechanism to explain. Complicated patterns, such as wide variety of profile configurations, rich in 
fluctuation or smooth structures, rapid rise vs. slower fall, weak precursor and secondary pulses 
etc, are common in various electrical disruptive discharges, e.g. various kinds of gas discharge, high 
vacuum disruption, or lightning flashes. McBreen et al. (1994) have noted that there is considerable 
similarity in the statistical properties of the light curves and peak flux distributions of GRBs and 
terrestrial lightning: The range in durations of the flashes and multiple stroke patterns within the 
flashes are well fitted by lognormal distributions and the integral number versus peak flux distributions 
by truncated lognormal distributions. Fishman et al. (1994a) discovered a dozen of brief, intense 
gamma-ray flashes of atmospheric origin, terrestrial gamma-ray flashes (TGFs), from the first two 
years observation of BATSE. Although TGFs have durations, typically a few milliseconds, much 
shorter than GRBs, their time profiles showing complex features in wide variety are quite similar with 
GRBs. TGFs have hard spectra with the hardness ratio H%2 = F3/F2 being about two times < H32 > 
for GRBs. For all bursts, of which the energy fluences F2 and F3 are available in the 3B catalog, we 
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Figure 15: Autocorrelation functions of light curves in four energy channels expected by the discharge 
model of Fig. 12 and Fig.14. The time lag is scaled by the time constant of discharging current decay. 

found < H32 >= 3.82. There exists an anti-correlation between GRB hardness and duration, hardness 
with short events in general being harder. Fig. 16 shows the relation between the hardness ratio H32 
and duration T90 in the 3B bursts and the regression line H32 = 5.1 — 1.1 lg(Tgo(s)). Extrapolating 
above relation to T90 = 5 x 10~ 3 s, we expect H32 = 7.6, just about two times < H32 >• Thus, besides 
the similarity in morphology, TGFs can be seen, in view of the spectral behaviour, as an extrapolation 
of GRBs. An evident correlation of TGFs with storm system indicates that these events are caused 
by electrical discharges to the stratosphere or ionosphere of the Earth (Fishman et al. 1994a), then 
TGFs give an observational evidence that the discharge process can produce high-energy bursts like 
GRBs. 

Another kind of high-energy transients similar with GRBs is the solar hard X-ray flare. It is not 
easy to make the distinction between GRBs and solar hard X-ray flares just by their time profiles. In 
burst recognition, solar flares distinguished from GRBs by their locations near the sun, their generally 
soft spectra, and their coincidence with solar X-ray events detected by the GOES satellite (Fishman 
et al. 1994b). Opposite TGFs, solar flares are in general longer and softer than GRBs. From 1991 
April 19 to 1995 November 28, there were a total of 4467 solar flares detected by BATSE, of which 677 
gave rise to a burst trigger (SDAC 1996). Yu et al. (1997) have found that the intensity distribution 
of BATSE solar flares is very similar with that of GRBs. Fig. 17 shows the intensity distribution of 
solar flares in the BATSE flare catalog and that of GRBs in the 3B catalog. The solar flare integral 
lgiV — lgP distribution also appears a —3/2 slope for bright events and saturated at low intensities 
like GRB does. After parallel translation the lg N — lg P diagram of solar flares is well coincidence 
with that of GRBs. This remarkable fact implies that the characteristic of GRB intensity distribution 
is not necessarily to be a measure of spatial homogeneity of burst sources, and that GRBs and solar 
flares may have similar production mechanism. Another noticed similarity between GRBs and solar 
flares is the delayed 1-GeV emission observed in both events (Kanbach et al. 1993; Hurley et al. 
1994). In fact, Alfven and Carlqvist (1967) suggested exploding discharges of electric double layers to 
be responsible for solar flares. 
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Figure 16: Correlation between hardness and duration of GRBs in 3B catalog. The arrow indicates 
the typical value of H32 for terrestrial 7-ray flashes detected by BATSE. 




Figure 17: Integral intensity distributions of GRBs and solar flares. Peak fluxes P of GRBs are 
measured on the 1.024s timescale for 3B bursts, that of solar flares are taken from the BATSE solar 
flare catalog in SDAC for the period from 1991 April 19 to 1995 November 28. The dashed line is 
drawn from the solar flare IgN — IgP distribution by a parallel translation lgiV + 0.25, lgP — 1.2. 
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